ON WEAKLY FORMULATED SYLVESTER EQUATIONS AND 

APPLICATIONS 



LUKA GRUBISIC AND KRESIMIR VESELIC 

Abstract. We use a "weakly formulated" Sylvester equation 

A 1 / 2 TM~ 1 / 2 - A~ 1 / 2 TM 1 ^ 2 = F 

to obtain new bounds for the rotation of spectral subspaces of a nonnegative 
selfadjoint operator in a Hilbert space. Our bound extends the known results 
of Davis and Kahan. Another application is a bound for the square root of 
a positive selfadjoint operator which extends the known rule: "The relative 
error in the square root is bounded by the one half of the relative error in 
the radicand" . Both bounds are illustrated on differential operators which are 
defined via quadratic forms. 



1. Preliminaries 

In this work we will study properties of nonnegative selfadjoint operators in a 
Hilbert space which are close in the sense of the inequality 



(1.1) \h{<t>,il>) - mfaVOI < vVmrnbP] 

where the sesquilinear forms h, m belong to the operators H, M respectively and 
m[ijj] = m(ip, V), h[4>] = h((f), (j>). 

In the first part of the paper we show that (|l.lfl implies an estimate of the 
separation between "matching" eigensubspaces of H and A. To be more precise 
one of the typical situations is: Let 

(1.2) < Ai(H) < A 2 (H) < • < A„(H) < D < A„ +1 (H) < ■ ■ • 

(1.3) < Ai(M) < A 2 (M) < • < A„(M) < D< A„ +1 (M) < • • 
be the eigenvalues of the operators H and M which satisfy ll.lfl then 



\\E H (D) E M (D)\\ < mm{ D _ MH) , ^a^M) I ^ 

Such an estimate 1 was implicit in 0. We then generalize this inequality to hold 
both for the operator norm || • || and the Hilbert-Schmidt norm ||| • \\hs- We also 
allow that En(D) and Em.(D) be possibly infinite dimensional 2 . 

In the second part of the paper we establish estimates for a perturbation of the 
square root of a positive operator. It will be shown that the inequality (|l.lfl implies 

\h2(4>,ip) - rri2{<i>,ip)\ < - v / /^p]m^J, 
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where the sesquilinear forms h,2,m2 belong to the operators H 1 / 2 ,]^ 1 / 2 , respec- 
tively. This will show that it is meaningful to consider weakly formulated Sylvester 
equations where all the coefficient operators are unbounded, cf. If .41 

Both of this problems will be solved through a study of the weak Sylvester 
equation, which reads formally 

(1.4) HI-XM = H 1/2 fM 1/2 . 

These two case studies represent two different classes of additional assumptions 
which have to be imposed on the coefficient operators H, M and F in order that 
(|1.4(l defines a meaningful operator X. 

The main novelty (and contribution) of this work is that we present an abstract 
study of the operator equation (|1.4(l in the case when only F is a bona fide operator. 
The expression H^i^M 1 / 2 need not possess an operator representation. In com- 
parison, was always a bounded operator for the Sylvester equations 
which were studied in 121 EH • 

Our first main result, contained in Theorem 12.11 below, extends our previous 
result from [7] in various ways. In particular, we allow the perturbed projection to 
be infinite dimensional. In the proof we also overcome a technical error contained in 
[Jj . We then extend this result to the case of other unitary invariant operator norms 3 
Particular attention is paid to the Hilbcrt-Schmidt norm because of its possible 
importance in applications. This special case is handled by another technique 
which allows an arbitrary interlacing of the involved spectra. 

The main object in this work shall be a closed nonnegative symmetric form in 
a Hilbert space. When dealing with symmetric forms in a Hilbert space, we shall 
follow the terminology of Kato, cf. jjj). For reader's convenience we now give 
definitions of some terms that will frequently be used, cf . E] ■ 

Definition 1.1. Let h be a positive definite form in 7i. A sesquilinear form a, 
which need not be closed, is said to be h-bounded, if Q(h) C Q{a) and there exists 
V > 

\a[u]\ < r)h[u] ueQ(h). 

If h is positive definite the space (Q(h), h) can be considered as a Hilbert space. 
The form a, which is /i-bounded, defines a bounded operator on the space (Q(h), h). 

Definition 1.2. A bounded operator A : Ti — > U is called degenerate if R(A) is 
finite dimensional. 

Definition 1.3. If H is a self adjoint operator and P a projection, to say that P 
commutes with H means that u € 2?(H) implies Pu e 25(H) and 

UPu^PUu, u6f(H). 

Definition 1.4. Let H and M be nonnegative operators. We define the order 
relation < between the nonnegative operators by saying that 

M < H 

if and only if Q(H) C Q(M) and 

\\M^ 2 u\\ < \\U^ 2 ul u e Q(H), 

3 Also called "cross-norms" in the terminology of |S| or "symmetric norms" in the terminology 

of nun]. 
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or equivalcntly 

m[u] < h[u], u € Q(h), 

when m and h are nonnegative forms denned by the operators M and H and 
M < H. 

A main principle we shall use to develop the perturbation theory will be the 
monotonicity of the spectrum with regard to the order relation between nonnegative 
operators. This principle can be expressed in many ways. The relevant results, 
which are scattered over the monographs 0IE], are summed up in the following 
theorem, see also [101 Corollary A.l]. 

Theorem 1.5. Let M = J A dEjvt(A) and H = j X dEn(X) be nonnegative oper- 
ators in TL and let M < H. Let the eigenvalues of H and M be as in iji.jj)) and 
ifO)) then 

(1) A e (M) < A e (H) 

(2) dim Eu("f) < dim Sm(7), for every 7 € K 

(3) A fe (M) < A fe (H), fc = l,2,--- . 

The infimum of the essential spectrum of some operator H is denoted by A e (H). 

With this theorem in hand we review spectral properties of operators H and M, 
for which there exists < e < 1 such that 

(1.5) (1 -e)m{u] < h\u] < (l + e)m[u], u G Q := Q{h) = Q(m). 
Let us assume h[u] > then m[u] > and 

(1.6) (1 —)h[u] < m[u] < (1 + v- —)h[u], 

1 — £ ' ' 1 — £ 

Inequality (|T3|) implies that N(H) = N(M), so (JTT3) holds for all u E Q. 

Lemma 1.6. Let m and h be nonnegative forms such that A e (M) > and A e (H) > 
and let $1.5p hold. Then 

(1.7) |Ai(H)-Ai(M)|<eAi(M) 

(1.8) |Ai(H) - A,(M)| < — Ai(H) 

1 — £ 

Ai(H) and Xi(M) are as in i)i.£J) and Assume that Aj_i(H) < Ai(H) < 

Ai+i(H) and 



(19) -^< m ax( A ' +l(H) - Al(H) A l( H)-A t _ l( H) y 

{ ' 1-e <maX lA. 1+1 (H) + A l (H)' A.^ + A^^H)' 1 / 

then 



(1.10) min = IMH)-yM)| < L 

a^m) A. 4 (H) Aj(H) 

//Ai_i(H) < A,(H) = ■ • ■ = A l+ „_!(H) < A 4+n (H) and 

1 ^ 1-e lA l+n (H) + A,(H)' A^ + A^H)' 1 / 
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then 

/- „, • |Aj-i(H)-A 7 -(M)| . , 

(1.12) argmin 1 \ M < i - 1 

(1.13) argmm' — 11 > 2 + n. 

jgn Ai+„(H) 

Proof. Estimates <|1.7ll - H1.8fl are a consequence of l|1.5|) - (|1.6|) and Theorem 11.51 
The rest of the theorem follows from a proof which analogous to the proof of pi 
Theorem 4.16]. We repeat the argument in this new setting. 
Let i j then 

|A,(H) - A,(M)| |Ai(H) - A,(H)| A*(H) + A,(H) |A,-(H) - \j (M) | A,(H) 



Aj(H) " Aj(H) + Aj-(H) Ai(H) A,(H) A*(H) 

A,-(H)\ s A,(H) 



- 7 V Xi(H)J 



> 



A,(H)7 l-EAi(H) 
|Ai(H)- Ai(M)| 



> 7 



A»(H) 

With this we have established (|1.1UI) . 1(1.12(1 - 1(1. 13(1 are a way to state (|1.1U|) in a 

presence of a multiple eigenvalue Aj(H). The proof follows by a repetition of the 
previous argument for j >i and j < i + n — 1. For instance, we establish (|1.12l) by 
proving 

lAj-iCHj-AjtM)! |Ai_i(H)- Ai_i(M)| 
Ai_i(H) > Ai-i(H) 
for all j > i. □ 

Remark 1.7. The significance of this Lemma is that it detects which spectral 
subspaces should be compared. When we were comparing discrete eigenvalues, the 
order relation between the real numbers (eigenvalues) solved this problem automat- 
ically. For spectral subspaces we need to assume more than ((1.5(1 in order to be 
able to construct meaningful estimates. Assumptions 1)1.9(1 and ((1.11(1 show how 
much more we (will) assume. 

Next we show that 1(1.5(1 implies 1(1. 1|) with 77 = e(l — e) -1 / 2 . To establish this 
claim we need a notion of a pseudo inverse of a closed operator. A definition from 
[P7| will be used. The pseudo inverse of a self adjoint operator H is the self adjoint 
operator Ht defined by 

V(ltf) = R(H)eX>(H) x , 

U^u + v) = H~y u e R(H), g e V(li) ± . 



It follows that Ht = H 1 in R(H). Note that we did not assume Ht to be bounded 
or densely defined. The operator Ht will be bounded if an only if R(H) is closed in 
Ti., see ^31- The operator Ht could have also been defined by the spectral calculus, 
since 

'0, A = 0, 

In (17j Weidmann has given a short survey of the properties of the pseudo inverse of 
a nondensely defined operator H. In particular, let Hi and H2 be two nonnegative 



Ht = /(H), /(A) = 
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operators in 2?(Hi) and T>(H.2) respectively then 

(1.14) ||H; /2 U || < HH^till * ||Hf Ml < ||H; /2t u ||. 

Analogously, let hi and hi be two closed, not necessarily densely defined, positive 
definite forms and let Hi and H2 be the self adjoint operators defined by hi and 
hi in Q(h±) and Q{hi). We say hi < hi when Q(h 2 ) C Q(/ii) and 

(1.15) hi[u] = ||H; /2 U || 2 < h 2 [u] = ||H 2 /2 u|| 2 , u G Q(M- 

Equivalently, we write Hi < H2 when hi < hi. Now, we can write the fact (11 . 14|) 
as 

(1.16) Hi<H 2 ^H 2 <Hj. 

In one point we will depart from the conventions in [Hj. 

Definition 1.8. A nonnegative form 

h(u,v) = (H^VH 1 / 2 ?;) 

will be called nonnegative definite when H^ is bounded. Analogously, the nonneg- 
ative operator H such that H* is bounded will also be called nonnegative definite. 

In the sequel we establish a connection between (|1.5f) and (|1.1|) when h and m 
are nonnegative definite forms. 

Lemma 1.9. Let H and M be nonnegative definite operators in a Hilbert space Ti. 
such that \1.5\) holds for < e < 1 . Let 



(1.17) 5 = H 1 / 2 M tl / 2 -Hti/2 M 1 /2 

then S is bounded and 

(i-i8) IW,^)| < -^=MH\\- 

VI -£ 

Proof. The closed graph theorem implies that the operator 



S = H^ 2 M n / 2 - HtV2MVa 

is bounded. Also, N(H) = N(M) = 1X1(5") and P N(S ) commutes with S. It is 
sufficient to prove the estimate for x, y G R(H). The assumption (|1.5fl gives 

\(h-m)(H^ 2 x,M^ 2 y)\ < e\\y\\ m^^x] 1 ' 2 . 

Analogously, 1)1.5(1 implies 

(1.19) || M l/ 2ti tl/2| < _L= , 



Altogether, the estimate 1(1.1811 follows. □ 

Now, we rewrite the conclusion of this lemma in the symmetric form setting. 
The result is given in the form of a proposition which we present without proof. 

Proposition 1.10. Let m and h be nonnegative definite forms and let there exist 
< e < 1 such that fO)) holds then N(H) = N(H) and 



\h(u, v) — m(u, v)\ < y/ h[u]m[v\. 

Vl - e 

When we only know that h and m satisfy p. 1(1 then we can establish a similar 
result about N(H) and N(M). 
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Proposition 1.11. Let m and h be nonnegative definite forms such that 11. lp holds 
then 

S = H^M* 1 / 2 - Hti/2MV2 
S* = Mti/2HV2 _ M^IT 1 / 2 

are bounded operators and \\S*\\ = \\S\\ < rj. Furthermore, N(H) = N(M) and a 
fortiori R(H) = R(M). 

The operator 5* has a special structure. Assume Mu = fiu and Hv — Xv, then 

(t>, Su) = A 1 ^ 2 (w, u)^ 2 - A~ 1 / 2 (u, u)^ 2 

(1.20) = A_^ (U;U ). 

vV 

The equation ll.2U[l suggests the distance function 

\X - m| 



which measures the distance between the eigenvalues of operators H and M. We 
state this result as the following corollary.. 

Corollary 1.12. Let Mit = fiu, \\u\\ = 1 and Hv = Xv, \\v\\ = 1 and let S be as in 
Proposition 1 1 . 1 il then 



Our theory is designed to be directly applicable to differential operators given in 
a weak form. This will enable us to obtain estimates for the difference between the 
spectral projections of the operators to which the theory of |2] does not apply, 
see Example 13 .41 below. 

2. Weak Sylvester equation 

Let us outline the general picture. We have an unbounded positive definite 
operator A and a bounded positive definite operator M. They are defined in, 
possibly, different subspaces of the environment Hilbert space H. Thus, Hm = 
R(M) is (of necessity) a closed subspace of TL and likewise 

D(AV2) n = R(A : / 2 ) = H A . 

Let the bounded operator F : Hm — > Ha be given, then we are looking for the 

bounded operator T : Hm — > Ha such that 

(2-1) 

(A^ 2 v, TM- 1 / 2 ^ - (A~ 1//2 v, TM 1 / 2 ^ = (v, Fu) , v£ V(A^ 2 ), u G Hm- 
Formally, we say that T solves the equation 
(2.2) AT - TM = A 1/2 FM 1/2 . 

Here G — A^^FM 1 ! 2 is naturally only a formal expression and does not represent 
a bona fide operator. In the case in which G be a bona fide operator equation (|2.2(l 
becomes the rigorous equation 

AT - TM = G, 
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called the (standard) Sylvester equation, cf . !B G] ■ The case when A and M are 
finite matrices has been considered in where (|2.2|l was called the structured 
Sylvester equation. 

We call the relation l|2.1|l the weak Sylvester equation. It represents a general- 
ization of the concept of the structured Sylvester equation 12.2fl from finite matrix 
setting to unbounded operator setting. The following theorem slightly generalizes 
the corresponding result from the joint paper jxj and corrects a technical glitch in 
one of the proofs. 

Theorem 2.1. Let A and M be positive definite operators in Tij^ and TLm, re- 
spectively and let F be a bounded operator from TIm into R(A X / 2 ) = Ha_. If M is 
bounded and 

( 2 - 3 ) W M W < p^J 

then the weakly formulated Sylvester equation 

(2.4) (A 1 / 2 *;, TM-^v) - (v, A^TM^u) = («, Fu) 

has a unique solution T, qiven by t(v,u) — (v,Tu) and 
1 f°° 

(2.5) t(v,u) = / {A 1 / 2 v, (A — i£ — d)~ 1 F(M — iC — d)^ 1 M 1 l 2 u)d( s , 

2tt J_ 00 

where d is any number satisfyinq 

(2.6) \\ M \\<d<-^—- . 



Proof. The uniqueness means that 

(2.7) (A 1 / 2 ?;, WM- l / 2 u) - (v, A^WM 1 / 2 ^ = 0, 

for u g %, v € P(A 1 / 2 ), has the only bounded solution W — 0. Let 

E n = / d E A i /2 (X), 
Jo 

then in particular 

(A 1/2 v,E n WM~ 1/2 v) - [y,A- 1/2 E n WM 1/2 uj = 0, 

for u G Hm , V G V(A 1 / 2 ) n E n H. Define the cut-off function 

. , s fx, D < x < n 
fn[x) = \ s 

In, n < x 

with D — 1/||A _1 ||. The operator /^.(A 1 / 2 ) is bicontinuous and 

(2.8) f n (A^ 2 )E n WM-^ 2 - f n {A 1 ' 2 )~ 1 E n WM 1 ' 2 = 0. 

Since /^A 1 / 2 ) and M 1 ^ 2 are bounded and positive definite operators, the standard 



Sylvester equation (|2.8|l has the unique solution 
(2.9) E n W = 0, neN. 

This is a consequence of the standard theory of the Sylvester equation with bounded 
coefficients, see The statement l|2.9[) implies W = 0. 
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Now for the existence. We use the spectral integral A = J A dE(X) to compute 

/OO pOO 
\\(A + i(-d)- 1 A 1 / 2 v\\ 2 d(= / {A 1 ' 2 v 1 \A-i(-d\~ 2 Av) d( 
-oo J — oo 

X d(E(X)A 1 / 2 v,A 1 / 2 v) 



oo poo 



d (\-d) 2 + C 2 



/ Xd(E(X)A 1 / 2 v,A 1 / 2 v) 
Jd J-, 

tt\ d{E{X)A 1 / 2 v,A 1 / 2 v) 



d( 



(X-d) 2 +( 

v, A 1 / 2 v) 

X-d 

(2.10) =Tt{A{A-d)- 1 v,v). 
Analogously, one establishes 

/oo 
|| (M d) -1 M 1/2 u|| 2 d( = Tt(M(d - M)~y u). 

-oo 

The convergence of these integrals justifies the following computation. Set 
1 f°° 

t(v, u) = / {A 1/2 v, (A - iC - d) _1 F(A/ - i£ - dy 1 M 1/2 u)dC, 

and then compute using (|2.10() and 12.1111 

1 r f°° i 2 

\r(v, u)\ 2 = _i_ [J ((A + iC - d)-^ 1 /^, F(M - iC - d)-^ 1 /^)^ 

|| (A + iC - A 1 / 2 ^!! || (M - i( - d^M^uWdC, 



< 11*11 



(2.12) < ^^(A(A~d)- 1 t>, W )(Af(d-M)- 1 u,w). 



(27T) 

4 

This in turn implies that the operator 

t(v, u) = (v, Tu) 

is a bounded operator and also gives the meaning to the formula (|2.5() . 
Now we will prove that this T satisfies the equation Note that 

A(A~ p-d)- 1 = I + (p + d)(A- p-d)-\ p£<r(A) 
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and then take v £ 2?(A) to compute 

i r 

~ ~2^ 



/OO 
(Aw, (A - iC - d)~ l F{M - if - d) _1 u) rfC 
-oo 

(v, (A - iC - d)~ 1 F{M - i( - d) _1 Mu) d( 



1 f v.p. / (v, F(M d)~ l u) d( 



2n 



« + d)((A - if - d) _1 », F(M - i( - d)- 1 !!) d( 
(i( + d)((A - i( - dy 1 v, F(M d( 

/"CO 

v.p. / ((A-iC-d^v^u) d( 



= (v,Fu). 

By a usual density argument we conclude that the operator T satisfies (|2.4H . □ 
Theorem 2.2. Let A, M and F be as in Theorem \2.1\ then 



\T\\ < 



D\\M\\ \\F\\ 



(D-d)(d-\\M\\) 2 
for any \\M\\ < d < D . The optimal d is d = (||M|| + D)/2 and then we obtain 



Proof. Estimate $FTty yields 



|r||<M|| A(A _ drl ||||M(d^)-' 11/1 ; ^ l|M|1 




(D-d)(d-\\M\\) 



This in turn implies the desired estimate. The optimality of the d = (\\M\\ + D)/2 
can now be checked by a direct computation. □ 

Remark 2.3. In fact, we will se that the estimate of Theorem l2.2l is optimal in the 
following sense. Let us consider the equation l|2.4|) in another light. Theorem 12.11 
gives a set of conditions when the equation (|2.1|l has a unique solution. Theorem 
12.21 then provides us with an estimate of this solution. 

Since for given F, under the conditions of Theorem 12. II there exists the unique 
T such that (|2.4() holds, we can define the so called "Sylvester operator" which 
associates the solution T to every operator F. The estimate H2.18II is then an 
estimate of the norm of the inverse of such an operator. 

The bound l|2.13|l is sharp in this sense as shows the following example. Let M 
and A be such that 

Mq=\\M\\q, Ap = Dp, 
for p and q one dimensional projections and let F = pq. 
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Then (|2.4(l is obviously satisfied by 

2.1. Allowing for a more general relation between er(M) and <r(A). An 
analogue of Theorem l2.1l holds. if the assumption 12.6(1 is replaced by a more general 
one, namely that the interval 

[IIM-i-MlMH] 

be contained in the resolvent set of the operator A. We omit the proof of the 
following result. 

Theorem 2.4. Let the operators A, M and F be as in Theorem \2. U and let their 
spectra be arranged so that 

(t(A) C (0,D_] U [D + ,oo), 
where < £>_ < || AJ- 1 1| - 1 ||M|| < D + . Then (in the sense of PT5)) ) 

T = / A 1/2 (A-i(-d)- 1 F(M-i(-d)- 1 M 1/2 dC 

27r J-oc 

+_L / a 1 /2( A - iC - g)- l F{M - K ~ gy'M^dt, 

27r J-oo 

where D- < g < ||M _1 || _1 and \\M\\ < d < D + , is the solution of the weak 
Sylvester equation j2.4\/ - We also have the estimate 

l|T|1 - l||M-i||-i-D_ + D + -\\M\\) m - 

2.2. Estimates in the Hilbert Schmidt norm. A bounded operator H : H — > 
7i is a Hilbert-Schmidt operator if fPiJ is trace class and then, cf. 8, Ch. X.1.3], 



(2.14) I H \\ HS := Trx/iFff. 

Let A and M be positive definite operators in Ha C H and 7Ym C H, re- 
spectively. We will analyze the weakly formulated Sylvester equation under the 
assumption that ||| F ||hs< oo and 

(2.15) gap( CT (M),a(H)) := inf > 0. 

AGct(A) 

To prove our result, we will need a basic result on the spectral representation of 
selfadjoint operators, see (181 Satz 8.17]. 

Theorem 2.5 (Spectral representation). For every selfadjoint operator H in a 
sparable Hilbert space 7i there exists a a -finite measure space (M,fi), a ^-measur- 
able function h : A4 — > R and a unitary operator V : 7i — > L 2 (A4,fi) such that 

H = V~ 1 H.V. 

Here H : L 2 (Ai , /i) — > L 2 (J\4 , /x) is the multiplication operator which is defined by 
the function h. 

We will also need the following theorem on the integral representation of Hilbert- 
Schmidt operators. For the proof see ^3 Satz 3.19]. 
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Theorem 2.6. A bounded operator T : L 2 (Mi,fi) — > L 2 {M.i,v) is a Hilbert- 
Schmidt operator if and only if there exists a function t G L 2 {M \ x M.i, /JXc) such 
that 

(Tg)(y) = / t(x,y)g(x)dfj, ^-almost everywhere, g£L 2 (M\,^). 

J Mx 

Furthermore, we have 

II T \\hS= \\A\l 2 (M 1 xM 2 4ixv)- 

We now prove a "Hilbert-Schmidt" version of Theorem l2.ll We will assume that 
II F \hs< 00 and that TL be separable. On the other hand, the spectra of A and 
M may be arbitrarily interlaced. 

Theorem 2.7. Let A and M be positive definite operators in TLa and TLm., re- 
spectively and let F : TLm. — * be a bounded operator. Assume further that 
II F \\hs< 00 an d gap(<r(M), c(H)) > then there exists a unique Hilbert-Schmidt 
operator T such that 

(2.16) (A 1/2 u,r]Vr 1/2 w) - (»,A- 1/2 TM 1/2 iij = (v,Fu) 

and 

(2J?) 1 T iHS - ga P (I(MWH))- 

Proof. The uniqueness of the bounded solution of the equation (|2.16|) follows by a 
double cut-off argument analogous to the one used in H2.8[) - l|2.9[) . We leave out the 
details. 

By Theorem 12.61 there exist measure spaces (Mm.,^) and (A4a,h), measurable 
functions m : .Mm — * K and a : Ma —> K and unitary operators U : TL — > 
L 2 {Mm,Ii) and V : TL — > L 2 (M a ,/j) such that 

A = V _1 A> 
M = C/^MZ/. 

Here we have taken A and M to be the multiplication operators which were defined 
by the functions a and m respectively. Since ||| F \\hs< 00 > the operator VFU : 
£ 2 (.MM,Ai) — * £ 2 (-Ma,m) is obviously a Hilbert-Schmidt operator and ||| VFU || 
|hs = 1 -P 1 |||j?s- We can therefore assume, without loosing generality, that we work 
with TLm = L 2 (Mm,v), 7~La = L 2 (Mm,v) and that A = A, M = M and F = 
VFU. 

Theorem 12 . 61 implies that there exists a function / S L 2 (Mm x Ma, (ixi/) such 
that 

(Fg)(y) = / f(x,y)g(x)dfj. i/-almost everywhere, g £ L 2 (Mm,p)- 

J Mrji 

Set 

fix y\ 

(2.18) t(a;,y) = afa)1/2 ' m(g;)1/2 , M x ^-almost everywhere. 

mix) 1 / 2 ~ a(y) 1 / 2 

Relation l|2.15|l and the positive definiteness of A and M imply that 

XOVM^ 1 

11 o(-)-m(-0 lli°°(^MxA4 A ,Mx v ) S gap((j(M)j(j(A)) 
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thus t G L 2 (Mm x Ma, M X v) and 

(2.19) ||i||L2 (A -! M x A -i AiA1 x 1 , ) < gap^dM) a (A)) H^C^mxjWa^x^- 
Now (|2.18[) can be rewritten as 

(2.20) a(y)^ 2 t(x, y)m(x)-^ 2 - a(y)-^ 2 t(x, y)m(x)^ 2 = f(x, y) 
The kernel t defines a Hilbert-Schmidt operator T with 



(v,Tu) — / v{y)t(x,y)u(x)dn dv. 



By taking integrals for v G V(A 1 ^ 2 ) and u G X^M 1 / 2 ) we establish that the equa- 
tion (|2.20j) is equivalent to (|2. 16j) and the estimate (|2.19jl implies (I2.17j) . □ 

2.3. Estimates by other unitary invariant operator norms. Let C(TL) be 
the algebra of all bounded operators on the Hilbert space H. We will consider 
symmetric norms ||| • ||| on a subspace S of C(7i). To say that the norm is symmetric 
on S C C{TL) means that, beside the usual properties of any norm, it additionally 
satisfies: 

(i) : If B e S, A, C G C{H) then ABC G S and 

I ABC |||< ||A|| I B I ||C||. 

(ii) : If A has rank 1 then ||| A \— \\A\\, where || • || always denotes the standard 
operator norm on C(7i). 

(iii) : If A G <S and U, V are unitary on H, then UAV G 5 and ||| [7j4F ||| = ||| 

(iv) : 5 is complete under the norm ||| • |||. 

The subspace S is defined as a ||| • |||-closure of the set of all degenerate operators in 
C{H). Such S is an ideal in the algebra C{H), cf. |HfTS|. Symmetric norms were 
used in [l] in the context of subspace estimates. If we assume, additionally to the 
assumptions of Theorem 13.21 that ||| F ||| < oo then there exists a unique bounded 
solution T of the weak Sylvester equation and 



JD\\M\\ 

~ D - \\M\\ 

We now prove this fact. 

Theorem 2.8. Let A and M be the selfadjoint operators which satisfy the assump- 
tions of Theorem 12.11 and let the symmetric norm ||| • ||| have the property 
(P) //sup I A n |||< oo and A — w-lim„A„ then A G S and 

II A |||< sup I A n If . 

V III F II < oo then there exists a unique bounded operator T such that 

(aV^TM- 1 /^ _ L,A- l ' 2 TM l l 2 u) = (v,Fu) 



m \/D\\M\\ ,„ 
T < — t; TT F 



L. Grubisic and Krcsimir Vcsclic 



13 



Proof. The proof follows by a cut-off argument. We (re)use the construction which 
was used in Q2.8J1 . Let f n (A x / 2 ) and E n be as in l|2.8[l . The equation 

(2.21) (.^(A 1 / 2 )^ T n M- 1/2 nj - (/^A 1 / 2 )-^, T n M^ 2 uj = (v, E n Fu) 
can now be written as the standard Sylvester equation 

f n (A^ 2 ) 2 T n - T n M = UA^EnFM 1 ^ 2 

which has the unique bounded solution T n : TLm ~ > R(E n ) and ||| T n |||< oo (this 
follows from [21 Theorem 5.2]). The operator E n T is bounded and satisfies the 
equation (|2.21|) therefore T n = E n T. Here we have tacitly assumed £(Ha) C C(TL). 
Furthermore, 

(2.22) A 1/2 E n TM- 1/2 - A- 1/2 E n TM 1/2 = E n F. 
We compute, using the property (i), 

II A-^EnTM 1 ' 2 II < 11 AJV211 II E n T II 

III A l / 2 E TM 1 / 2 III > V_S^lHE! Ill e T III 

From these estimates and Q2.22t we obtain the uniform upper bound 
JD\\M\\ \/D\\M\\ 

Since E n T — > T in the strong operator topology, Property (V) and the uniform 
bound (|2.23|l imply ||| T |||< oo and the desired estimate follows. □ 



3. Perturbations of spectral subspaces 

When comparing two spectral subspaces of operators H and M, which satisfy 
(|1.1|) . we have to make an additional assumption on the location of the spectra. 
Namely we assume that there exist D\ < Di such that the interval [Dx,D-i\ C K 
is contained in the resolvent sets of both H and M. Let Q = En(Di) and P = 
Em{D\). We want to estimate the norm of P — Q. The following description of a 
relation between a pair of orthogonal projections in a Hilbert space will be sufficient 
for our considerations. For the proof see 

Theorem 3.1 (Kato). Let P and Q be two orthogonal projections such that 

||P(I-Q)||<1. 

Then we have the following alternative. Either 

(1) R(-P) and R(Q) are isomorphic and 

||P(I-Q)|| = ||Q(I-P)|| = ||P-Q|| or 

(2) R(P) is isomorphic to true subspace of R(Q) and 

||Q(I-P)|| = ||P-Q||=1. 
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To ease the presentation set P± = I — P and Q± = I — Q. First, let us consider 
the case when h and m are positive definite. With the help of Proposition II . 1 II we 
shall later reduce the nonnegative definite case to the positive definite one. 

We define the operators 

(3.1) A = Q ± iJQ ± , H = QHQ, M = PMP and W = P±MP±. 

We shall not notationally distinguish the operators A, M, W and H from their 
restrictions to the complement of their respective kernels. Obviously, 

H = + A, M = Af + W 

and we compute, for S from (|1.17(l . 

Q ± SP = {H 1/2 Q ± PM- 1/2 - B.- 1/2 Q ± PM 1/2 )P 

= A 1/2 Q ± PA'r 1/2 - A- 1/2 Q ± PM 1/2 

(3.2) = A 1/2 TM- 1/2 - A- 1/2 TM 1/2 . 

Here we have defined T = Q±P. If we assume that dim(Q) = dim(P) < oo then 
Theorem 13.11 yields 

\\P-Q\\ = \\T\\. 

The case when dim(Q) = dim(P) = oo will follow in a similar fashion. 
The operator equation can be written in the following variational form 

(3.3) (A 1/2 v, TM-^u) - (A-^ 2 v, TM x ' 2 u) = (v, Su), 

v e V(A 1/2 ), u e R(P), 
which we have called the weakly formulated Sylvester equation. 

Theorem 3.2. Let the positive definite forms m and h be given such that 
holds. Let there exist D\ < D2 such that the interval [Di,D2] C R be contained in 
the resolvent sets of both H and M. Set Q = En(Di), P = Em(-Di) and assume 
n < {D 2 -D 1 ){D 2 D 1 )- 1 I 2 then 



(3-4) l|P-QII<-^^. 

Proof. T = Q±P is the unique solution of the equation (|3.3|) . Theorem 12 . 21 implies 

IITII < 



D 2 Xn(M) 



\j (D 2 -d)(d-\ n (M))- 

for any A n (M) < d < D 2 . The optimal d equals (D2+A 2 " (M)) and since ||A/|| < D 1 

we conclude 

y/£> 2 A n (M) VZhLX 

~*~ g^A^(M) Th^Di 

Analogous argumentation for T = P±Q, with the roles of H and M in H3.3(l being 
interchanged, yields the inequality 

II^QII<^#^<i. 

L>2 — U\ 

Theorem ()3.1|l now implies that 

||QxP|| = ||PxQ|| = ||Q-P||. 
This in turn establishes 1)3. 4|l . □ 
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In the case in which h is only nonnegative definite, assumption implies that 
N(M) = N(H) and R(H) = R(M), since H and M are selfadjoint. This in turn 
allows us to conclude that TV := R(P) n N(H) c R(Q), so 

Q = Q-P^, P = P-Ph 
are orthogonal projections and 

||Q-P|| = ||Q-P||. 

Since R(P) C R(H) and R(Q) C R(H) we can reduce the problem to the positive 
definite case. 

Theorem 3.3. Let the positive definite forms m and h be given such that 
holds. Let there exist < L\ < L2 < D\ < D2 such that the intervals \L\, L2] C K 
and [Di,Z?2] C K be contained in the resolvent sets of both H and M. Set Q = 
E-a [Li , L2] , P — -Em[-Di, D2] and assume 



y/Lhl\ VL2L1 



D 2 -D 1 



Li 



r/ < 1 



then 



(3.5) 



l^-QII < 



VD2D1 , jL^Lx 



D2-D1 



L2 — L\ 



V 



Proof. The assumption L\ > implies that we may assume, without losing any 
generality that we have the positive definite forms m and h. Theorem 12 . 71 and the 
same argument as in Theorem 13 . 21 implies 



\\p±Q\\ = \\Q±p\\ 

This in turn allows us to conclude that 



\p-QW 



(3.6) 



l^-QII < 



D 2 -Dt L 2 - L x 



□ 



Numerical experiments with the Sturm-Liouville eigenvalue problem, which were 
performed in |7j, illustrated that in some situations the results of Theorems 13.21 
and 13.31 yield considerably sharper estimates of the perturbations of the spectral 
subspaces than the results of We now show that our theorems also apply in 

situations in which the theory from does not. 

Example 3.4. Take H, M as selfadjoint realizations of the differential operators 

d d d d 

q(x) — , /3(x) — , 

dx dx ' dx dx' 

respectively, in the Hilbert space Ti — L 2 (I), L a (finite or infinite) interval with, 
say, Dirichlct boundary conditions and non-negative bounded measurable functions 
a(x), (3(x) which satisfy 



Now, the form 



\(3{x) — a(x)\ < rj\f (3{x)a(x) 



S(u, v) = h(u, v) — m(u, v) 
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is not — in general — representable by a bounded operator. This rules out an appli- 
cation of the subspace perturbation theorems from [UEj. On the other hand both 
of our Theorems 13.21 and 13.31 apply and yield the corresponding estimates, e.g. 



L>2 — L)\ 

when we know that [D\, D2] is contained in the resolvent sets of both H and M. 

Theorem 12 . 71 can also be applied to yield a Hilbert-Schmidt version of Theorems 
IO and 1331 

Theorem 3.5. Let the positive definite forms m and h be given such that 1 1.1}) 
holds. Assume P and Q are projections which commute with the operators H and 
M respectively and let A, M, W, H as in fO)) . If ||| Q±SP || H s< 00, ||| P ± S*Q \\ 
|hs< 00 ond 

gap«A), <r(M)), gap(a(W),a(H)), 
are both positive. Then Q±P, P±Q and P — Q are Hilbert-Schmidt operators and 

(3-7) III Q,P Ills < j^rj^ 



(3-8) I P±Q § 



3 p(a(A),a(M)) 2 

QSP± Ills 



lllls - ga P (<7(W),(7(ff))2 

(3-9) lP-Qlls< 



2 ^ ill QxSP Ihs , I QSP X \\ls 



p(a(A),a(M)) 2 gap(<r(W), a(H)) 2 ' 

Proof, by construction (|3.1() the operator T = T\ = Q± satisfies the Sylvester 
equation l|3.3|l . which in this setting has the form, cf. I|2.16|l . 

(3.10) (A 1 / 2 v, TM-Wu) - (A^^v, TM^ 2 u) = (v, Q±SPu), 

veV{A 1/2 ) 7 ueV(M 1/2 ). 

On the other hand, the operator T = T2 = P±Q satisfies the "dual" equation, cf. 
Proposition 11.111 

(3.11) (W^ 2 v, TH-^ 2 u) - (W-^ 2 Vl TH^ 2 u) = {v, P ± S*Qu), 

veV(W 1/2 ), ueV{H 1/2 ). 

Now, 

(P-Q) 2 = Q±P + P±Q, 
where by Theorem 12 1 71 both Q±P and P±Q are Hilbert-Schmidt 4 and 

\P-Q lis = MQxP + P±Q) = Tr(PQj_P) + Tr(QPj_Q) 

=1 Q±p lis + III Px-Q lis 

Using PTT7|) . we see that estimates (|3~?|) - (|3~§j) hold. □ 

Corollary 3.6. Let the positive definite forms m and h be given such that \1.1}) 
holds. //II S |||hS< 00 an d the other conditions of Theorem \3. 5\ hold. Then 

III <? Ill 2 

a 12) hip -qui 2 < 111 J 111 hs 

{ - ' 111 VlllHS - min{gap( C r(A),a(M)) 2 ,gap(a(W),a(F)) 2 } ' 



4 To prove this equality one can use the singular value analysis from Alternatively, one 
could use the property (V) from Theoreml2.i 
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Proof. Just note that 

III Q±sp f HS + i QSP± ll 2 H s<|5|^ s . 

□ 



4. Further properties of the operator S — an application in the 

NUMERICAL ANALYSIS 

We will now present an application of Theorem 12.21 in numerical analysis. This 
will also demonstrate a role played by the new Hilbert-Schmidt norm estimates. 

Assume now that we are given a positive definite operator H such that (|1.2|) 
holds. Let P be an orthogonal projection such that R(P) C Q(h) and dim R(P) = n. 
We aim to obtain estimates of 

(4-1) I E n (D)- P||| 

for I • |||= || • || and ||| • ||| = ||| • \\ H s- 

We estimate l|4.1f) by an application of Theorcm l2.2l feauivalentlv Theorem 13. 5fl . 
Theorem l2 . 21 will allow us to improve 7 Theorem 3.2] inasmuch as that we establish 
estimates for the Hilbert-Schmidt norm and not just the spectral norm. 

The properties of the main perturbation construction from [7] , cf . [21 E] , will be 
summarized for reader's convenience. 

We start by defining the positive definite form 

(4.2) h P {u, v) = h(Pu, Pv) + h(P ± u, P ± v) 

and a self adjoint operator Hp which represents the form hp in the sense of Kato. 
It was shown, see [5] El E] that 

(1) the form hp is positive definite , hence there exists the positive definite 
operator Hp which represents hp in the sense of Kato. 

(2) Q(h) = Q(h P ) 

(3) R(P) reduces H P . 

(4) H" 1 — Hp 1 is a degenerate self adjoint operator. 

(5) Let Sh p :—h~hp and let SH P be the bounded self adjoint operator such 
that 

(4.3) («, SHfv) - Sh p (H p 1/2 u, H p 1/2 v) 

then SH P is a degenerate operator and d\mR(SH P ) = 2n. 

(6) The values 

(aa\ ■ f (^H-V)-(^HpV) X 

(4.4) ^ = max mm \ 

dimS—n—i 

together with their negatives are all non-zero eigenvalues of SH P . Further- 
more, rji are all the singular values of the operator 8Hf P. 

(7) 

(4.5) \Sh p ((f>, VOI < VnVhrW^W] 

The estimates from |7| Theorem 3.2] only use information which is contained in 
■q n . New theory allows us to take advantage of other 7^. 
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Proposition 4.1. Let P and hp be as in \4- 6 4) and let- 



then 



S = H 1/,2 Hp 1//2 - H-V2 H ^ /2 



SQ II < SH * Q ■ 



S III : 



V 1 - Vn 

III SHf II 



\/l - Vn 

Here SHf is the degenerate operator from \4-S\j , Q is any projection, rj n is given 
by \4-4\l and III ' III * s any unitary invariant norm. 

Proof. 

(*£,SQ<£) = tt p (H- 1 / 2 ^,H p 1/2 Q0) = ^ p (H p 1/2 (H^ /2 H- 1 / 2 )VAH p 1/2 g0) 

= (V,, (H^U-^YSHf Q<j>), 

(|4.5|) and ljl. 19|) imply ||Hp 2 H~ 1/2 | < 1/^/1 - ??„. Property (i) of the symmetric 
norm ||| • ||| and the fact that 8Hf Q is a degenerate operator allow us to complete 
the proof. □ 

This proposition leads to an improved version of (7J Theorem 3.2]. Observe that 
I 8H^ I depends only on r\i from l|4.4|l . 

Theorem 4.2. Let h be as in and let P and hp be as in \4-2jj and T)i as in 

Q. Set 

Dp := max 



^GR(P) \\ljj\\ 2 

and assume rj n (l — r/n) -1 < (D — Dp)(D + Dp) then 

(4. 6) , EK( > n)PL ,< i^m . 

D - Dp y/l-rin 
Here ||| • ||| is any unitary invariant norm which has Property (P). 

Proof. Set T = (En(A n ))xP and apply Theorem 12. 81 to estimate the norm 

I (En(K))±P I- PropositionlOnow implies g3, cf. Corollary EHfl 2, Corollary 

3.1] and [3 Proposition 6.1]. □ 

Assume ||| • ||| = ||| • \\hSi then Theorem 14.21 yields the estimate 



(A7\ i (w (\ w phi < ^ DDp vM±IE±M 

( 4 -< j III (-kH(A„)) ±J P jj S < — . 

D - Dp \J\-r\n 

Remark 4.3. Note that ||(Sh(A„))±P|| = \\E H (\ n ) - P\\, cf. Theorem O A 
similar relation holds for a general unitary invariant norm since according to 
Corollary 3.1] and |2 Section 2] we have ||| (Eu{\ n )) ± P ||| = ||| Pj_E H (X n ) ||| and 

(4.8) I E u {X n ) - P 1=1 (£ H (A n ))xP + P±^h(A„) ||| . 

An estimate of (|4.8f) is obtained by a combination of Proposition ^. II and available 
(depending on an application) information on the separation of the involved spectra, 
cf. Corollary 13.61 We have not specified a general estimate on ||| En(X n ) — P ||| 
since we consider such estimates to be highly application dependent and we would 
not like to prejudice their form. 
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N 


5 


6 7 8 9 10 


1 (-Eh(A 2 ))±P \\hs 


4.4e-3 


2.0e-3 l.le-3 6.0e-4 3.7e-4 2.4e-4 


\/hDp N y/rfi + r% 
A3 ~ D Pn y/\ - r/2 


2.2e-2 


1.0e-2 5.3e-3 3.3e-3 2.2e-3 1.5e-3 


Vs 1 (R$) + s 2 (R») 
A3 — Dp N 


2.0e-2 


1.4e-2 9.6e-3 7.2e-3 5.5e-3 4.4e-3 



Table 1. Error estimate from Theorem 14 . 21 and the true error 



We will now evaluate 1|4.7[) on the example from Section 4]. There we have 
considered the positive definite operator H which is defined by the symmetric form 

h(u, v) = (u'v' — auv) dt 







u, v e {/:/,/' E L 2 [0, 2tt], e i8 /(0) = /(2tt)} = V(h). 
The eigenvalues and eigenvectors of the operator H are 



(±*+^) -a, z± k (t)=e- i ( ±k+ £) t , ken 



2 



In standard notation we have 

Ai(H)=wo, A 2 (H) = w_ 1) A3(H)=wi, 

ui = z , u 2 = z-i, u 3 = zi. 

For numerical experiments we chose 8 = ir — 10~ 4 and a = 0.2499 so that the 
eigenvalues Ai and A2 are "small" and tightly clustered. As a test space we chose 
yfj — spanju;^, W2}, where and u>^ are generated by the smooth N point 
equidistant cubic interpolation of the known eigenfunctions u± and u 2 . Take Pn 
such that R(P N ) = y%. Since y% C D(H) both Theorem [O] and the bounds from 
P] apply. Set = H(j)+ [<f>, H<p)4>, for \\<fr\\ = 1. Since , € 2?(H) we conclude 
that r w N and r w N are bona fide vectors. Set 

~{r w N ,r w N ) {r w N,r w N) 
The competing bound from 2 is 



(4.9) I (E Ui x 2 )) ± p N i\ hs < vEm^m . 

A3 - D Pn 

On Table Q we have displayed the actual measured error in the first line, in the 
second line we display the bound from 1)4. 7JI and in the third line Davis-Kahan 
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bound (|4.9|l . We see that with the improvement of the approximation the advantage 
of the bound from Theorem 14.21 over l|4.9|) grows, see Tabled Further examples, 
where a numerical advantage of l|4.6|l over (|4.9|) is more stunning, are given in [7|- 
We repeat the results of the numerical experiments from [7] on Table There we 
try to estimate the approximation error in the vector in the || • ||-norm by an 
application of Theorem 13.21 Otherwise the makeup of Table El is the same as the 
makeup of Tabled 



N 


6 7 8 9 10 


\\E u {Xi) - Pn\\ 


2.0e-3 l.le-3 6.0e-4 3.7e-4 2.4e-4 


^2 - dp N ^/l - T] 2 


1.5e0 6.2e-l 3.5e-l 2.2e-l 1.5e-l 


A2 — dp N 


3.6e+2 2.1e+2 1.5e+2 l.lc+2 8.9c+l 



Table 2. Approximations for u\ (here we use dp N :— min^ e R(p N ) -rr^w) 



We now present a variation on this example where (|4.9(l does not apply whereas 
(|4.6|) still gives useful information. We chose = spanjZ^, l^}, where and 
I2 are generated by the N point equidistant continuous linear interpolation of u\ 
and U2 then r t N and are no longer bona fide vectors. Subsequently, l|4.9|) does 
not apply any more but Theorem 14.21 is still applicable. Take now Qn such that 
R(Qjv) = yh- The results are presented on Table 

The performance of the bound l|4.6|l is influenced by the quotient 

I^p w -A 2 | 
Dp N 

Dp N is an approximation 5 of A2 and in this example we have measured 
\D P -A9I 

■ L -^ ^>0.17, N = 100, 120, 140. 

°Pn 

The (under)performance of the bound (|4.6(l correctly detects this approximation 
feature of R(Q N ), cf. Table 

5. Estimates for perturbations of the square root of a nonnegative 

operator 

In this section we will show that there are interesting applications of the equation 
(|2.1(l even when all of the coefficients A, M and F are unbounded. To demonstrate 



To be more precise Dp N is Rayleigh— Ritz approximation to A2(H) from the subspace R(Pjv)- 
For more on the Rayleigh— Ritz eigenvalue approximations see 0. 
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N 


100 120 140 


1 {E n (X 2 ))±Q N \\hs 


5.2024e-005 3.6126e-005 2.6541e-5 




8.7374e-003 6.9293e-003 5.7302e-003 





Table 3. 



this we will generalize the known scalar inequality 6 

(5.1) j= — < , h,m>0. 

ymh lymh 

to positive definite Hermitian matrices or, more generally, to positive, possibly 
unbounded, operators in an arbitrary Hilbert space. One of the obtained bounds 
is 

(5.2) ||Af-V 4 (M-V a - ff- 1 / 2 )!/- 1 / 4 !! < i||M- x / a (M - 2f)if -1 / 2 ||. 
In ^2] a related bound for finite matrices was obtained. It reads 

(5.3) \\H- l l\M- 1 ' 2 - H-^H- 1 ^ < | + 0( V 2 ), 

This is a more common type of estimate — the error is measured by the "unper- 
turbed operator" only — while in our estimate the error is measured by H and M 
in a symmetric way. The latter type of estimate is convenient, if both operators H 
and M are known equally well and we are interested in a possibly sharp bound. Our 
bound is obviously as sharp as its scalar pendant. It is also rigorous, in contrast 
to (|5.3|) which is only asymptotic. Moreover, l|5.2|l will retain its validity for fairly 
general positive selfadjoint operators in a Hilbert space. The bound (|5.2I) is a "rel- 
ative bound" which may be convenient in computing or measuring environments 
(cf. related bounds obtained for the eigenvalues and eigenvectors of the Hermitian 
matrices in and the literature cited there) . Also, this bound is readily expressed 
in terms of quadratic forms, which will be convenient for application with elliptic 
differential operators as will be shown below. 

The idea of the proof is very simple, especially in the finite dimensional case 
which we present first, also in order to accommodate readers not interested in 
infinite dimension technicalities. 

The basis of our proof is the obvious Sylvester equation (cf. |15p 

(5.4) M 1 ' 2 ^ 1 ' 2 - H 1 ' 2 ) + (A/ 1 / 2 - H 1 ' 2 )!! 1 ' 2 — M — H, 

valid for any Hermitian, positive definite matrices H and M. We rewrite this 
equation in the equivalent form 

(5.5) M X I A XH- X I A + M-^Xff 1 / 4 = T 



"The relative error in the square root is bounded by the half relative error in the radicand" . 
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with 

(5.6) T = M- l ' 2 {M-H)H- 1 ' 2 , X = M^M" 1 / 2 - H- 1 / 2 )^ 1 ^, 
which is immediately verified. This equation has a unique solution 

poo 

(5.7) X = / e'^^M-'^TH-^e-^^dt. 

Jo 

(just premultiply (|5.5|l by e~ M 1 2 *M _1 / 4 , postmultiply by e~ H 1,2t H~ 1 l i , inte- 
grate from to oo and perform partial integration on its left hand side). Hence for 
arbitrary vectors <p, ip we have 

\{X^^)\ 2 < \\T\\ 2 (^J™ ^e-^'^M-V^M^^H^^ndtj 

P OO OO 

< l|r|| 2 / We-^^M'^UW^t / We-^'^H-^^fdt 
Jo Jo 

(5-8) = ^IIIHlW, 



where we have used the obvious identity 

l Cdt = 

2 



(5.9) / e~ 2Ct Cdt = \l 



n 



for C = H- 1 / 2 , M- 1 ' 2 , Thus, (J52J holds true. 

We now turn to the Hilbert space Ti. of arbitrary dimension. We assume that H 
and M are positive selfadjoint operators. This implies that all fractional powers of 
H and M are also positive. Neither of these operators need be bounded (or have 
bounded inverse). 

Theorem 5.1. Let H and M be positive selfadjoint operators in a Hilbert space 
X having the following property (A): I?(H 1/ ' 2 ) = X^M 1 / 2 ) and the norms 
and HM 1 / 2 • || are topologically equivalent. Then the same property is shared by 
H 1 / 2 andM 1 / 2 . The operators 



M -i/2 H i/ 2 MV2H-V2, M-V4H1/4, MV^H- 1 /*, 



(5.10) H" 1 /2M 1 /2 ; H 1 /2M- 1 /2 ; H- 1 /4 M 1 /4 j h^m- 1 ^ 

are well defined and bounded. Let 



(5.11) T = M 1 /2 H - 1 /2_ M - 1 /2H 1 /2 

and 



(5.12) X = MV4H-1/4 - M-^h 1 / 4 
then 

(5.13) 11*11 < ^irni. 
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Proof. The fact that the square roots inherit the property (A) is a consequence of 
Lowner type theorems (see e.g. |Sj, Ch.V, Th. 4.12). The corresponding pairs of 
operators in (|5.1U|I are mutually adjoint e.g. M -1 / 2 !! 1 / 2 = H 1 / 2 ]^ -1 / 2 etc. Obvi- 
ously, l|5.11|l and i|5.12[l reduce to T, X from l|5.11|l . if the space is finite dimensional. 
The equation l|5.5[l becomes here 

(5.14) (IH- 1/4 !i, M 1/4 v) + (IH 1/4 «, M" 1/4 «) = {Tu, v) 

for ueV A = ^(H 1 / 4 ) n ^(H- 1 / 4 ) and similarly for v and M. We will now prove 
this. 

The left hand side of (|5.14(1 equals 

(M^h-iMh- 1 /^ M 1 / 4 ^) - (H~ 1,/4 w, HV^M^TTJm 1 / 4 ^) 



+ (H 1 / 4 u, H- 1 /^ 1 /^- 1 /^) - (M _1 / 4 H 1 /4H 1//4 ti, M.~ 1 l i v) = 
(U-^ 2 u,M 1/2 v) - («,«) + («,«) - (H 1 / 2 !*, M _1 / 2 u) = 
(M 1 / 2 H _1,/2 m, u) - (M- 1 / 2 H 1 / 2 m,w) = (Tu,«). 
Now, substitute in l|5.14|l 

(5.15) v = e- M " 1/at M-VV, u = e - H " 1/Jt H- 1 ^ 

for any G X^M" 1 / 2 ), V € ^(H^ 1 / 2 ). Note that 

M-V^DJM-V 2 ) and H _1//4 I?(H _1 / 2 ) 

are invariant under e~ M / *, e _H 1 *, respectively so, in (|5.15|) we have u G Pa 
and u € 2?m- Then integrate (|5.15(l and use partial integration: 

f {Xe-^'^n- 1 ' 2 ^ e- M_I/2 V)* + f (Xe- H_I/2 V, e^^'M" 1 / 2 ^ = 
Jo Jo 

_ f(il e - H " /J V, e- M " /2t 0)^ + f Ve- H " 1/a V, e- M " /2 *M- 1 / 2 0)^ = 
Jo a* Jo 



+ /Ve- H " /2 ^, e - M ~ 1/2t M- 1 / 2 0)^ = 
Jo 

^,0)-(X e - H - 1/2 >,e- M - 1/2 ^) = 

/ S (Te- H " /2t H- 1 / 4 ^,e- M " /2 *M- 1 / 4 ( /))dt. 
Jo 

In the limit s — > oo by using the functional calculus for H, M, respectively and 
monotone convergence for spectral integrals we obtain 

in the norm. Hence 

POO 

(5.16) (XV, 4>)= (Te- H_1/2t H- 1 /V,e- M " 1/2t M- 1 / 4 0)^ 

Jo 

where the integral on the right hand side is, in fact, Lebesgue as shows the chain of 
inequalities in (|5.8|) which are valid in this general case as well. Here the identity 
(15.91) is used in the weak sense: 

(e- 2Ct C<j>, 4>)dt = ((/), (f>)/2, <j> G V(C) 
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for any positive selfadjoint C. Thus, 

|(^,0)| 2 <||T|| 2 (^»(0»/4. 

□ 

Remark 5.2. The main assertion (|5.13|) of Theorem l5.1l is obviously equivalent to 
the following statement: the inequality 



\m(cf),ip) - h(4>,i>)\ < e^h(4>, <j>)m{4>, ^) 

implies 



|m 2 (^, ip) - hi{(j>, if>)\ < -VMS (j))m2(ip,tp) 

where the sesquilinear forms ft., to, /i2,to,2 belong to the operators H, M, H 1 / 2 , 
M 1 / 2 , respectively. Thus, our theorem will be directly applicable to differential 
operators given in weak form. 

Example 5.3. Let again H and M be as in Example 13.41 That is to say take H, 
M as selfadjoint realizations of the differential operators 

d d d d 

dx dx ' dx dx 

in the Hilbert space TL = L 2 (I) (again I can be a finite or infinite interval) with 
the Dirichlet boundary conditions and non-negative bounded measurable functions 
a(x), (3{x) which satisfy 



\P{x) — a(x)\ < e^(3(x)a(x) 

Now 

|(M 1/2 0,M 1/2 ^)-(H 1/2 0,H 1/2 ^)| 2 < Qf \p(x) -a(x)\\ij;'tx)<t>'(x)\dx^ < 

e 2 a{x)W{x)\ 2 dt J^f3{x)\(j>'{x)\ 2 dt = e 2 ^ 1 ' 2 ^ 2 ^ 1 ' 2 ^ 2 
hence ||T|| < e and Theorem 15.11 applies yielding 

or, equivalently, in the terms as in Remark 15. 21 



\m 2 ((f>,ip) - h 2 ((j) 1 ip)\ < -^Jh 2 {(j),4>)m2(^-,^)- 
6. Conclusion 

With this work we complete our study of the weak Sylvester equation which 
started in 0. A notion of a weak Sylvester equation was introduced in [7] as a tool 
on a way to obtain invariant subspace estimates for unbounded positive definite 
operators. With this paper we show that there are applications of the concept of 
a weak Sylvester equation outside the theory of Rayleigh-Ritz spectral approxi- 
mations. We have extended out theory to infinite dimensional invariant subspaces 
and have obtained estimates of the difference between the corresponding spectral 
projections in all unitary invariant norms. With this results we have developed a 
counterpart of the sin O theorems from pQ for perturbations of operators which are 
only defined as quadratic forms. 
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Due to the very singular nature of integral representations (which can not be 
avoided by reformulation of the integrals) of the solution to the equation H2.4JI . cf. 
formula (|2.5[) . we were not able to extend the technique from £Q to prove that in 
the setting of Theorem 12.71 assumption ||| F ||| < oo also implies that there exists a 
bounded solution T such that ||| T |||< oo. We believe that this statement is true, 
but the proof will have to remain a task for the future and will most likely require 
another technique. 

An application of the concept to a perturbation of the square root of a positive 
definite operator shows that there are other application areas for weakly formulated 
operator equations and that the developed techniques are (and hopefully will be) 
easily adaptable to new situations. The applications which we have reported in this 
paper are presented as an illustration only. Further applications will be the subject 
of the future work, cf. 
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